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ABSTRACT 

Tangled, primordial cosmic magnetic fields create small rotational velocity perturba- 
tions on the last scattering surface (LSS) of the cosmic microwave background radi- 
ation (CMBR). Such perturbations can contribute significantly to the CMBR tem- 
perature and polarization anisotropies at large I > 1000 or so, like the excess power 
detected by the CBI experiment. The magnetic contribution can be distinguished from 
most conventional signals, as they lead to CMBR polarization dominated by the odd 
parity, B-type signal. Experiments like DASI and WMAP have detected evidence for 
CMBR polarization at small I. Many experiments will also probe the large I regime. 
We therefore calculate the polarization signals due to primordial magnetic fields, for 
different spectra and different cosmological parameters. A scale-invariant spectrum of 
tangled fields which redshifts to a present value B — 3 x 10 -9 Gauss, produces B-type 
polarization anisotropies of ~ 0.3 — 0A/J.K between / ~ 1000 — 5000. Larger signals 
result if the spectral index of magnetic tangles is steeper, n > —3. The peak of the 
signal shifts to larger / for a lambda-dominated universe, or if the baryon density is 
larger. The signal will also have non-Gaussian statistics. We also predict the much 
smaller E-type polarization, and T-E cross correlations for these models. 

Key words: magnetic fields-cosmic microwave background-cosmology:theory-large- 
scale structure of Universe. 



1 INTRODUCTION 

The origin of large-scale cosmic magnetic fields remains an 
intriguing question. It is quite likely that magnetic fields in 
astronomical objects, like galaxies, grew by turbulent dy- 
namo action on small seed magnetic fields (cf. Ruzmaikin, 
Shukurov & Sokoloff 1988; Beck et al 1996). However, this 
idea is not without difficulties in view of the constraints im- 
plied by helicity conservation and the more rapid growth 
of small-scale magnetic fields (Cattaneo & Vainshtein 1991; 
Kulsrud & Anderson 1992; Gruzinov & Diamond 1994; Sub- 
ramanian 1998, 1999; Blackman & Field 2000; Kleoorin et 
al 2000; Brandenburg 2001; Brandenburg & Subramanian 
2000; Brandenburg, Dobler & Subramanian 2002). Magnetic 
fields with larger coherence scales may also be present in 
clusters of galaxies (Clarke, Kronberg & Bohringer 2001) 
and at high redshifts (Oren & Wolfe 1995). Such large- 
scale coherent fields could present further problems for the 
dynamo paradigm. Alternatively, galactic or cluster fields 
could largely be a remnant of a primordial cosmological 
magnetic field (cf. Kulsrud 1990), although, as yet, there 
is no entirely compelling mechanism for producing the re- 
quired field. They could be possibly generated during a pe- 



riod of inflation or at a phase transition, perhaps with an 
almost scale-invariant spectrum (Turner & Widrow 1988; 
Ratra 1992; cf. Grasso & Rubenstein 2001 for a review). 
A primordial field, whose present-day strength is of order 
10 -9 Gauss, and is tangled on galactic scales, can also af- 
fect the process of galaxy formation (Rees & Reinhardt 1972; 
Wasserman 1978; Kim, Olinto & Rosner 1996; Subramanian 
& Barrow 1998a, SB98a hereafter). It is of considerable in- 
terest, therefore, to find different ways of limiting or detect- 
ing such primordial fields (see Kronberg 1994 and Grasso & 
Rubenstein 2001, Widrow 2003 for reviews). 

The nearly isotropic nature of the CMBR already places 
a limit of several nano Gauss on the present strength of any 
uniform (spatially homogeneous) component of the mag- 
netic field (Barrow, Ferreira & Silk 1997). Observations of 
CMBR anisotropies also provide potentially powerful con- 
straints on tangled magnetic fields (Subramanian & Barrow 
1998b (SB98b), 2002 (SB02)). Such fields produce vortical 
perturbations, which are overdamped in the radiation era 
and can then survive Silk damping (Silk 1968) on scales 
much smaller than the compressional modes (Jedamzik, 
Katalinic & Olinto 1998; SB98a). So their signal, if present, 
will be particularly evident at small angular scales below the 
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conventional Silk damping scale or at multipoles I > 1000 or 
so (SB98b, SB02). Intriguingly, the CBI experiment identi- 
fied significant excess power in the CMBR anisotropy spec- 
trum up to I = 3500 (with 2 sigma limits of 14 — 31 fjK at 
/ > 2010) (Mason et al 2002). We argued in SB02 that tan- 
gled magnetic fields which redshift to a present-day value 
of -Bo = 3 x 10~ 9 Gauss (see below for a definition of Bo), 
can contribute a non-negligible fraction of this signal. Other 
possibilities include the Sunyaev-Zeldovich effect (Bond et al 
2002; Komatsu & Seljak 2002), primordial voids (cf. Griffiths 
et al 2002) or features in the primordial power spectrum (cf. 
Cooray & Melchiorri 2002). One needs to isolate and test 
for the possible magnetic field contribution. We point out 
here that this can be done by looking for the corresponding 
distinctive polarization signals. 

Indeed, Seshadri and Subramanian (2001; Paper I) 
showed that tangled magnetic fields create distinctive small- 
scale (large I) polarization anisotropy, dominated by the odd 
parity B-type signals. Mack et aZ.(2002; MKK02) gave es- 
timates for such signals, but in the low / < 500 regime. 
Recently the DASI experiment detected E-type polarization 
and both DASI and WMAP have detected the T-E cross 
correlation, albeit at smaller I ~ 300 (Kovac et al 2002; 
Kogut et al 2003). Several experiments to probe the polar- 
ization anisotropy at the large I regime are also underway or 
being planned (VSA, ACBAR, ATCA, Polatron and Planck 
Surveyor) . Motivated by these experiments, we compute the 
B-type polarization signals at larger values of I, for a wider 
variety of cosmological parameters and spectral indices, than 
were made in Paper I. In the next section, we first recapitu- 
late the arguments of Paper I. We then present in Section 3, 
approximate analytic estimates of the signals. The B-type 
anisotropy for various models, obtained from a detailed nu- 
merical integration is given in Section 4. We also calculate 
there numerically the E-type polarization anisotropy and the 
T-E cross correlations for these models. Our predictions will 
allow comparison with future observations and help in de- 
tecting or ruling out significant magnetic field contributions 
to the signal on small angular scales. 



2 THE POLARIZATION ANISOTROPY 

Polarization of the CMBR arises from the Thomson scat- 
tering of radiation from free electrons, and is sourced by 
the quadrupole component of the CMBR anisotropy. The 
evolution equations for the moments, Oi, Ei and Bi, of the 
temperature anisotropy (AT/T), the electric (E-) type and 
the odd parity, magnetic (B-) type polarization anisotropics, 
respectively, for vector perturbations, have been derived in 
detail by Hu & White (1997a; HW97) (see also Paper I and 
SB98b). For vector perturbations, the B-type contribution 
dominates the polarization anisotropy (HW97). We there- 
fore give details of its calculations and summarize the re- 
sults for the E-type contribution and the T-E cross corre- 
lations. The quadrupole anisotropy source term for polar- 
ization is given by P(k, r) = [O2 — v6-E2]/10. Here k is the 
co-moving wave number, r the conformal time. One can ana- 
lytically estimate P using the tight-coupling approximation, 
kLj(r) <C 1, where Z/ 7 (r) is the co-moving, photon mean 
free path. First, to leading order in this approximation, we 
have zero quadrupoles, and a dipole 61 = vb, where «s(fc, r) 



is the magnitude of the rotational component of the fluid ve- 
locity vf , in Fourier space. However, to the next order the 
quadrupole is not zero. It is generated from the dipole at the 
'last but one' scattering of the CMBR. From Eq. (60), (63) 
and (64) of HW97, we get 6 2 = -AEh/VE = 4fcL 7 « s /(3 v / 3) 
and hence P = 9 2 /4 = kL y v B / '(3^3). Using this in Eq. (77) 
and (56) of HW97 gives an estimate for Bi and the angular 
power spectra Cf B due to B-type polarization anisotropy 
(cf. Paper I), 
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Here ji(z) is the spherical Bessel function of order I, and 
to the present value of r. The 'visibility function', g(ro,r), 
determines the probability that a photon reaches us at epoch 
tq if it was last scattered at the epoch r. We adopt a flat 
universe throughout, with a total matter density fi m and a 
non-zero cosmological constant density = 1 — £l m today. 
We now briefly recall the arguments detailed in Paper I. 

Firstly, we approximate the visibility function as a 
Gaussian: g(r ,r) = (27r<7 2 ) -1/2 exp[-(r - t,) 2 / (2a 2 )], 
where r» is the conformal epoch of "last scattering" and 
a measures the width of the LSS. To estimate these, we 
use the WMAP results (cf. Spergal et al 2003), that the 
redshift of LSS z, = 1089 and its thickness (FWHM) 
Az = 194. To convert redshift into conformal time we use, 
r = 6000ft -1 ((a + a^) 1 / 2 — alq 2 )/Qm 2 , valid for a flat uni- 
verse (cf. Hu & White 1997b). Here, the expansion fac- 
tor a = (1 + z)' 1 and a eq = 4.17 x 10~ 5 (f2 m ft 2 ) -1 (ft is 
the Hubble constant in units of 100 km s -1 Mpc -1 ). For 
the A-dominated model suggested by WMAP results, with 
n m = 0.27, Qa = 0.73, we get r» = 201. 4ft -1 Mpc and 
a — 11.5ft -1 Mpc. For an Q, m = 1 model, with the same 
baryon density, (cf. tl b = 0.0224ft" 2 ), we use the expressions 
given in Hu and Sugiyama (1995) to estimate r* = 131.0ft -1 
Mpc, and a = 8.3ft -1 Mpc. We use these numbers in the nu- 
merical estimates below. 

To evaluate Cf B , one needs to estimate vb, the rota- 
tional velocity induced by magnetic inhomogeneities. The 
magnetic field is assumed to be initially a Gaussian random 
field. On galactic scales and above, the induced velocity is 
generally so small that it does not lead to any appreciable 
distortion of the initial field (Jedamzik, Katalinic & Olinto 
1998, SB98a). So, the magnetic field simply redshifts away 
as B(x, t) — bo(x)/a 2 . The Lorentz force associated with 
the tangled field is then Fl = (V x bo) x bo/(47ra 5 ), which 
pushes the fluid and creates rotational velocity perturba- 
tions. These can be estimated as in SB02 or Paper I, by 
using the Euler equation for the baryons. On scales larger 
than the photon mean-free-path at decoupling, where the 
viscous effect due to photons can be treated in the diffusion 
approximation, this reads (SB02) 
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Here, p 7 is the photon density, p b the baryon density, and 
77 = (4/15)p 7 Z 7 the shear viscosity coefficient associated 
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with the damping due to photons, whose mean-free-path 
is Z 7 = (n e CTT) _1 = L y a(t), where n e is the electron den- 
sity and ar the Thomson cross-section. We have also ig- 
nored here a metric perturbation term which is subdom- 
inant at large I (cf. Paper I). For z» ~ 1089, we get 
L 7 (-r«) ~ 1.83/" 1 Mpc, where f b = (fi b /t 2 /0.0224) (the 
WMAP value). We have defined the Fourier transforms 
of the magnetic field, by bo(x) = ]T] k b(k) exp(ik.x) and 
F(k) = E p [b(k + p).b*(p)]p-[k.b*(p)]b(k + p). The pro- 
jection tensor, Pij(k) = [8ij — kikj/k 2 ] projects F onto its 
transverse components perpendicular to k. 

The comoving Silk damping scale at recombination, 
Ls = kg 1 ~ 10 Mpc, separates scales on which the radiative 
viscosity is important (kLs 1) from those on which it is 
negligible (kLs <C 1). For kL s <ti 1, the damping due to the 
photon viscosity can be neglected compared to the Lorentz 
force. Assuming negligible initial rotational velocity pertur- 
bation, we can integrate the baryon Euler equation to get 
vf = GiD, where d = 3P ij F j /[16np ] and D = r/(l + S») 
(see SB02 or paper I). Here po is the present-day value of 
p 7 , and S, = (3p b /4p 7 )(r») ~ 0.59/f,. For kL s » 1, we 
use the terminal-velocity approximation, neglecting the in- 
ertial terms in the Euler equation, to balance the Lorentz 
force by friction. This gives vf = Gi(k)D, but with now 
D — (5/k 2 Lj), on scales where diffusion damping operates. 
The transition Silk scale can also be estimated by equating 
vf in the two cases, to give k.s ~ [5(1 + 5»)/(rL 7 (r))] 1 ^ 2 . 

To compute the Cf B s we need to specify the spec- 
trum of the tangled magnetic field, say M(k). We define, 

< 6i(k)6j(q) >= 8 K(i Pij(k)M(k), where <5 kjq is the Kro- 
necker delta which is non-zero only for k = q. This gives 

< b 2 , >= 2f(dk/k)Al(k), where A 2 b (k) = k z M(k)/(2-K 2 ) 
is the power per logarithmic interval in k space residing 
in magnetic tangles, and we replace the summation over k 
space by an integration. The ensemble average < \vb\ 2 >, 
and hence the Cf B s, can be computed in terms of the mag- 
netic spectrum M(k). It is convenient to define a dimension- 
less spectrum, ra(k) = A 2 (k)/(B 1 /2), where Bo is a fiducial 
constant magnetic field. The Alfven velocity, Va, for this 
fiducial field is, 
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where B-g = (£>o/10~ 9 Gauss). We will also consider as in 
SB02, power-law magnetic spectra, M(k) = Ak n cut-off at 
k = fe c , where k c is the Alfven- wave damping length-scale 
(Jedamzik, Katalinic & Olinto, SB98a). We fix A by de- 
manding that the smoothed field strength over a "galactic" 
scale, fee = 1/iMpc -1 , (using a sharp fe-space filter) is Bo, 
giving a dimensionless spectrum for n > — 3 of 

m(fc) = (n + 3)(fc/fc G ) 3+ ". (4) 



Here, A 2 = k' A < |t> B (fc, r*)| 2 > /(27r 2 ) is the power 
per unit logarithmic interval of k, residing in the net ro- 
tational velocity perturbation, and R, — to — r*. In the 
opposite limit, ka > 1, we get 1(1 + l)Cf B /(2n) ^ 
(kL 1 (r*)/3) 2 (^/4:)(Al(k,T,)/(ka)\ k=l/Rt . At small wave- 
lengths, Cf B is suppressed by a 1 /ka factor due to the finite 
thickness of the LSS. Further in both cases, the polariza- 
tion anisotropy, ATf B (l) w (fcL 7 (r»)/3) x AT (I), where, 
AT(l) is the temperature anisotropy computed in SB02. 
We can now put together the above results to derive ap- 
proximate analytic estimates for the CMBR polarization 
anisotropy induced by tangled magnetic fields. As a mea- 
sure of the anisotropy we define the quantity ATf B (l) = 
[1(1 + l)Cf B /2tt] 1/2 T , where T = 2.728 K is the CMBR 
temperature. On large scales, such that kL s < 1 and ka < 1, 
the resulting CMBR anisotropy is (see Paper I) 

AT# s (0 = To(-) 1/2 I(k) k2L ^ Tt)VlTt 
32 
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Here, I = kR* and we adopted the A-dominated model fa- 
vored by WMAP, with n A = 0.73, Vl m = 0.27, Q. b h 2 = 
0.0224 and h = 0.71 (in Paper I, we used a purely matter- 
dominated £l m = 1 model). We also use the fit given by Hu 
& White (1997b) to calculate r = 6000fe _1 ((l + a eq ) 1/2 - 



alq 2 )(l — 0.0841 ln(fi m ))/f2m^, valid for flat universe. On 
scales where kLs > 1 and ka > 1, but fcL 7 (r») < 1, we get 
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The function I 2 (k) in the Eqs.(5) and (6) is a dimensionless 
mode-coupling integral given in our previous papers (cf. Eq. 
(7) of Paper I). Analytic approximations to I(k) exist for 
power-law spectra and for k <§C k c (as generally relevant 
even at high Z; cf. Paper I and MKK02) For n > -3/2, 



t2/u\ - 28 (n + 3) 2 . k 3 k c 3+2n 
' W ~ 15(3 + 2n) ( fe ) { k^> ' 



(7) 



dominated by the cut-off scale k c . For n < —3/2 (SB02), 

I 2 (k) = 8 -(n + 3)(^r 2n (8) 

independent of k c , where we neglect a subdominant term 
of order (k c /kf +2n < 1; note that k < k c . A nearly 
scale-invariant spectrum, say with n = —2.9, then gives 
ATf B (l) ~ O.WpK(l/1000) 2 1 for scales larger than the 
Silk scale, and ATf B (l) ~ 0.51pK(l/2000y 0A , for scales 
smaller than Ls but larger than L 7 . Larger signals result 
for steeper spectra, n > —2.9 at the higher I end. 



3 ANALYTIC ESTIMATES 

The dominant contributions to the integral over r in Eq. (1) 
come from a range a around the epoch r — r* . Furthermore, 
ji(k(ro — t)) picks out (k, r) values in the integrand which 
have k(ro— t) ~ I. Thus, following the arguments detailed in 
Paper I and SB02, for to < 1 we get the analytical estimate, 
1(1 + l)Cf B /(2tv) « (fci 7 (r*)/3) 2 (^/4)A 2 (fc,r»)U =i/flt . 



4 NUMERICAL RESULTS 

We now compute ATf B (l) for the above spectra, by eval- 
uating the t and k integrals in Eq.(l) numerically. For this 
we first express Eq. (1) explicitely in terms of the magnetic 
correlation function. We have for I > 1, (ATf B (l)/T ) 2 = 
J™(dk/k)(l 2 kVZl(k)U(k)/V8) 2 , where 
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Figure 1. AT0 B versus I predictions for different cosmological 
models and magnetic power spectrum M(k) cc fc n , for _B„g = 3. 
The bold solid line is for a standard flat, A-dominatcd model, 
with n A = 0.73, n m = 0.27, Q b h 2 = 0.0224, h = 0.71 and al- 
most scale invariant spectrum n = —2.9. The long dashed curve 
obtains when one changes to n = —2.5, while the short dashed 
curve is for a larger baryon density Q),h = 0.03. The dotted curve 
gives results for a Q m = 1 and Qjy = model, with n = —2.9. 
These curves show the build up of power in B-type polarization, 
due to vortical perturbations from tangled magnetic fields which 
survive Silk damping at high I ~ 1000 — 3500. The eventual flat- 
tening or slow decline is due to the damping by photon viscosity, 
although this is only a mild effect as the magnetically sourced 
vortical mode is overdamped. We also show for qualitative com- 
parison (dashed-dotted curve), the B-type polarization anisotropy 
due to gravitational lensing, in the canonical A-CDM model, com- 
puted using CMBFAST (Seljak & Zaldarriaga 1996; Zaldarriaga 
& Seljak 1998). The signal due to magnetic tangles dominate for 
I larger than about 1000. Finally, the thin solid line gives the ex- 
pected galactic foreground contribution estimated by Prunet et 
al (1998), which is also smaller than the predicted signals. 
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In doing the above integral numerically, we retain the ana- 
lytic approximations to I(k) and D with a transition be- 
tween the limiting forms of D, at wavenumber ks- (The 
function ji is treated as in codes like CMBFAST). The nu- 
merically evaluated results are shown in Figures 1 and 2. 

We see that for Bo ~ 3x 1CF 9 G, this leads to a predicted 
RMS B-type polarization anisotropy in the CMBR of order 
0.3 — 0.4/iif for 1000 < I < 5000, for a nearly scale-invariant 
power law spectra with n = —2.9. Larger signals result, at 
the high I end, even for a moderately steeper spectral index, 
with n — —2.5, which has more power on small-scales (com- 
pare the long dashed and solid curves). Much larger signals 
result from even steeper spectra with n > —2.5. However, 
these cases are probably ruled out for B-g ~ 3, as they 
lead to an over production of gravitational waves (Caprini 
& Durrer 2002). Hence we do not display the results for 



such spectra explicitely. The peak of the polarization signal 
shifts to larger values of I with increasing A (compare the 
solid and dotted curves), due to an increase in R*. This also 
happens for larger baryon density (compare the solid and 
short dashed curves), due to the decrease in L 7 and hence 
the damping effects of radiative viscosity. The /-dependence 
of ATp B (l) got by analytic approximation matches very well 
with that obtained by numerical integration. The amplitude 
is, however, somewhat overestimated by the analytic treat- 
ment; by a factor of about 1.2 and 1.6 respectively, for the 
regimes above and below the Silk scale. So one can indeed 
get a reasonable idea of the signals from the analytics, but 
for more accurate amplitudes, one needs the numerical inte- 
gration done above. 

We have also numerically computed the E-type polar- 
ization (C EE ) as well as the T-E cross correlation {Cf E ) 
for the above models. For this we use Eqns. (77), (79) and 
(56) given in HW97 and evaluate the k and r integrals nu- 
merically, adopting the same analytic approximations for 
I(k) and D given above. In calculating the T-E signal, we 
have included a small polarization contribution to T as well, 
which contributes negligibly to T, but significantly to Cf E . 
The cross correlation due to the vortical modes induced by 
magnetic tangles has a negative sign and so we define the 
corresponding effective 'temperature' after taking a modu- 
lus of Cj E . We show in Figure 2 the resulting E and T-E 
anisotropics, for the standard A-dominated model. For com- 
parison, we have also give the temperature (T) and B-type 
polarization anisotropics (for T we have corrected a normal- 
isation error made in SB02). The E-type polarization has a 
peak value of ~ O.l/i K and so is much smaller that the B- 
type signal, as expected for vector perturbations (cf. HW97). 
The T-E cross correlation is ~ 0.1 - 0.2/j.K, for I > 1000. 
However, both E and T-E power are subdominant to that 
produced by the standard scalar perturbations. 



5 DISCUSSION 

We have re-examined the small angular scale polarization 
anisotropy induced by tangled magnetic fields, for different 
cosmological parameters and spectral indices. A major mo- 
tivation arises from ongoing and future experiments, which 
probe this large I regime. Tangled magnetic fields generate 
vortical perturbations, which lead to a distinctive polariza- 
tion anisotropy dominated by a B-type contribution. Vor- 
tical perturbations survive Silk damping on much smaller 
scales than do compressional modes and their damping due 
to the finite thickness of the LSS is also milder. By con- 
trast, in the standard non- magnetic models the polarization 
anisotropy is dominated by E-type contributions. A scale- 
invariant spectrum of tangled fields which redshifts to a 
present value of Bo = 3 x 10~ 9 Gauss, produces B-type 
polarization anisotropies of order ~ 0.3 — OAfiK between 
I ~ 1000 — 5000. Larger signals are produced for steeper 
spectra with n > —3. These signals are also larger than the 
small B-type polarization induced by gravitational lensing 
at I > 1000 or the expected galactic foreground contribution. 
The peak of the signal shifts to larger !, in a A-dominated 
universe, or in a universe with larger O b . Our results comple- 
ment the small I (I < 500), and purely analytical estimates 
of MKK02. The polarization anisotropy peaks or troughs 
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Figure 2. The predicted anisotropy in temperature (dotted line), 
B-typc polarization (solid line), E-type polarization (short dashed 
line) and T-E cross correlation (long dashed line) up to large 
I ~ 5000 for the standard A-CDM model, due to magnetic tangles 
with a nearly scale invariant spectrum. 

could be much larger, because the non-linear dependence 
of Cf B on M(k), implies non-Gaussian statistics for the 
anisotropics. Clearly, with the sub-micro-Kclvin sensitivities 
expected from experiments like the Planck Surveyor, these 
signals can be detected. 

An important contributor to small scale anisotropies, 
like that seen the CBI experiment, would be the Sunyaev- 
Zeldovich effect. This signal can be isolated by its frequency 
dependence. Also, scattering in clusters produces a much 
smaller statistical polarization anisotropy compared to the 
magnetically-induced signals (even individual clusters pro- 
duce maximum signals of only ~ 0.1/iK (cf. Sazonov & 
Sunyaev, 2000)). The polarization induced by the SZ ef- 
fect, primordial voids, or features in the power spectrum 
are all expected to be predominantly E-type, in contrast 
to predominantly B-type signals predicted here. Primordial 
magnetic fields can also lead to depolarization due to differ- 
ential Faraday rotation. This effect is only important at fre- 
quencies lower than about 16.4GHz(B_9/3) 1 ^ 2 (Kosowsky 
& Loeb 1996; Harari, Hayward & Zaldariaga 1997). Addi- 
tional B-type polarization can arise from tensor modes gen- 
erated during inflation or those induced by the primordial 
magnetic fields (cf. Durrer, Ferreira & Kahniashvili 2000; 
MKK02; SB02); but these are expected to be important 
only at / < 100 or so. Helicity in the magnetic spectrum can 
also leave interesting signatures on the CMBR (Pogosian, 
Vachaspati & Winitzki 2002). 

In summary, sensitive observational searches for B-type 
polarization anisotropies at large I will allow us to detect or 
constrain primordial, tangled magnetic fields. They will also 
tell us whether such fields are significant contributors to the 
excess power at large I detected by the CBI experiment and 
help probe possible new physics in the early universe. 
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